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Abstract. This review article gives an overview of recent progress in the field of 
non-equilibrium phase transitions into absorbing states with long-range interactions. 
It focuses on two possible types of long-range interactions. The first one is to replace 
nearest-neighbor couplings by unrestricted Levy flights with a power-law distribution 
P(r) ~ r - d - a controlled by an exponent a. Similarly, the temporal evolution can 
be modified by introducing waiting times At between subsequent moves which are 
distributed algebraically as P(At) ~ (At) -1 ~ K . It turns out that such systems with 
Levy-distributed long-range interactions still exhibit a continuous phase transition 
with critical exponents varying continuously with a and/or k in certain ranges of the 
parameter space. In a field-theoretical framework such algebraically distributed long- 
range interactions can be accounted for by replacing the differential operators V 2 and 
dt with fractional derivatives V 7 and 9 t K . As another possibility, one may introduce 
algebraically decaying long-range interactions which cannot exceed the actual distance 
to the nearest particle. Such interactions are motivated by studies of non-equilibrium 
growth processes and may be interpreted as Levy flights cut off at the actual distance 
to the nearest particle. In the continuum limit such truncated Levy flights can be 
described to leading order by terms involving fractional powers of the density field 
while the differential operators remain short-ranged. 
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1. Introduction 



In non-equilibrium statistical physics the study of phase transitions from a fluctuating 
active phase into absorbing configurations has been a very active field in recent 
years [1]. The primary goal of these studies is to classify continuous phase transitions 
by their universal properties. The most important universality class of absorbing 
phase transitions, which plays a paradigmatic role as the Ising model in equilibrium 
statistical mechanics, is Directed Percolation (DP) [2-5]. This class of phase transitions 
is characterized by three independent critical exponents f3, u±, and u\\ which depend 
solely on the dimensionality of the system but not on the microscopic details of the 
dynamics. Directed percolation corresponds to a specific type of field theory which first 
appeared in the context of hadronic interactions [6-8]. The structure of this field theory 
led Janssen and Grassberger to the so-called DP-conjecture, stating that any model that 
complies with a certain set of conditions should belong to the DP class [9, 10]. One of 
these conditions is that the interactions have to be local in space and time. 

The present review focuses on the question how the universal properties of 
absorbing phase transitions change if one generalized these models as to include non-local 
interactions. There are several motivations to study such models. On the one hand, 
directed percolation is often used as a simple model for the spreading of epidemics 
among a population of spatially distributed individuals. In order to make the modelling 
more realistic, various generalizations have been introduced involving infections over 
long distances, incubation times, as well as other memory effects [11-15]. On the other 
hand, DP-like processes with long-range interactions seem to play a role in certain 
models for non-equilibrium wetting, where an interface is grown on top of an inert 
substrate according to dynamic rules that violate detailed balance. Finally, from a 
more fundamental point of view it is of theoretical interest how long-range interactions 
can be accounted for in a field-theoretic description of non-equilibrium phase transitions. 

A possible generalization of directed percolation with a non-local transport 
mechanism was proposed already 40 years ago by Mollison in the context of epidemic 
spreading [12]. Motivated by empirical observations he suggested to consider a model 
where the disease is transported in random directions over long distances r which are 
distributed according to a power law 

P(r)~r- d ~ a , (1) 

where a > is a control exponent. In the literature such algebraically distributed long- 
range displacements are known as Levy flights [16] and have been studied extensively 
in the literature, e.g. in the context of anomalous diffusion [17-20]. It should be 
emphasized that a does not introduce an additional length scale, rather it changes 
the scaling properties of the underlying transport mechanism. It turns out that Levy- 
distributed interactions in models for epidemic spreading do not destroy the transition, 
instead they change the critical behavior, provided that a is small enough [21-23]. 
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More recently, Janssen et al showed that such transitions with spatial Levy flights 
can be described in terms of a renormalizable field theory [24], which allows one to 
compute the critical exponents to one-loop order. Moreover, this field theory predicts 
an additional scaling relation so that only two of the three critical exponents turn out to 
be independent. These results were confirmed numerically by Monte Carlo simulations 
in Ref. [25]. 

In the continuum limit Levy flights can be generated by certain non-local linear 
operators called fractional derivatives. Although fractional calculus is an interesting 
subject on its own (see e.g. [26-29]) this review is not intended to discuss mathematical 
issues in this direction, instead it is meant as a practical guide for statistical physicists 
explaining how these concepts can be used in the context of absorbing phase transitions. 

Likewise it is possible to introduce a similar long-range mechanism in temporal 
direction. Such 'temporal' Levy flights can be understood in terms of algebraically 
distributed waiting times At that may be interpreted in the context of epidemic 
spreading as incubation times between catching and passing on the infection. The 
probability distribution of these waiting times is given by 

P(At) ~ Ar 1- " , (2) 

where k > is the temporal Levy exponent. Unlike spatial Levy flights, which are 
isotropic in space, such temporal Levy flights are directed forward in time in order 
to respect causality. Recently a field-theoretic renormalization group calculation of 
DP with temporal Levy was presented in Ref. [30], while the mixed case of epidemic 
spreading by spatial Levy flights combined with algebraically distributed incubation 
times was considered in Ref. [31]. As a further step towards more realistic models, 
Brockmann and Geisel [32] recently studied epidemic spreading by Levy flights on a 
inhomogeneous support in the supercritical phase. 

Another type of long-range interactions is motivated by models for non-equilibrium 
wetting, where an interface grows on top of a hard-core substrate [33-39]. In these 
models it is particularly interesting to study the effective dynamics of pinning sites 
at the transition from a bound to an unbound state. It turns out that in some 
regions of the parameter space the pinning sites can be interpreted as active sites of 
a directed percolation process. However, the effective rate for binding a neighboring 
site (corresponding to local infections in the language of DP) may depend on the actual 
configuration of the adjacent detached segment of the interface, giving rise to a varying 
rate that depends algebraically on the distance to the nearest active site. In other 
regions of the parameter space, where the interface fluctuates so strongly such that it 
may reattach far away from other pinning sites, non-local infections with power-law 
characteristics are observed. However, these non-local interactions differ from Levy 
flights in so far as they are truncated at the actual distance of the nearest pinning site. 
This truncation leads to strikingly different features, as will be discussed in SectJU 
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The present review is divided into three parts. The following section briefly reviews 
some fact about long-ranged random walks. Sect. [3] discusses the absorbing phase 
transitions with unrestricted Levy flights in space and algebraically distributed waiting 
times, while truncated Levy flights are considered in Sect. HI 



2. Levy flights, fractional derivatives, and anomalous diffusion 



Before considering interacting particle systems such as directed percolation with long- 
range interactions, let us briefly recall some facts about long-ranged random walks 
by Levy flights, also known as anomalous diffusion [40]. We will consider two different 
types of long-range effects, namely, isotropic Levy flights in space as well as algebraically 
distributed waiting times. 



2.1. Spatial Levy flights 

A Levy flight in d dimensions is a random displacement r — > f + s according to a 
probability distribution which for large distances decays algebraically as 

P(s) ~ |s|- d ~ CT . (3) 

Here the exponent a is a control parameter which determines the asymptotic power-law 
characteristics of long-range flights. Note that this parameter does not induce a length 
scale, instead it changes the scaling properties without breaking scale invariance. The 
distribution P(s) has to be normalized by J d d sP(s) = 1 which implies that a > 0. 
Moreover, the normalization requires a lower cutoff of the power law ([3]) in the limit of 
small s, for example in form of a minimal flight distance. 

Let us first consider a particle which moves by a series of spontaneous uncorrelated 
Levy flights. The statistics of such a long-ranged random walk is described by the 
probability density p(f, t) to find this particle at position r at time t, which may also 
be interpreted as the particle density in an ensemble of many non-interacting particles 
performing Levy flights^- Obviously, p(r, t) evolves in time according to the integro- 
differential equation 

J^(f, t) oc J d d s P(s) [p(f + s, t) - p(f, *)] . (4) 

In order to write this equation in a more compact form, let us introduce a linear non-local 
operator V CT which acts on a function f(r) by 

V CT m = / \s\- d - [f{?+ s) - f(fj\ , (5) 

I As we will see in subsection 2.2 the interpretation of p(r, t) as a particle density is no longer valid in 
systems with algebraically distributed waiting times. 
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where 



7T 



i/2 r(-f; 



is a cr-dependent normalization constant. Neglecting higher-order derivatives caused by 
the short-range cutoff, the equation of motion can be written as 

J^(f,t) = V>(f,t) (7) 

which has the same structure as the ordinary diffusion equation apart from the fact that 
the Laplacian has been replaced by V CT . 

In the literature the operator V CT is known as a fractional derivative because it has 
certain algebraic properties that generalize those of ordinary derivatives. For example, 
it is straight forward to show that for < a < 2 the action of the operator V°" on a 
plane wave amounts to bringing down a prefactor of the form 

V a e* f = -\k\*e**. (8) 

Therefore, in momentum space the equation of motion ^p(k,t) = — \k\ a p(k, t) implies 
that each mode with wave number k decays exponentially as e~' fc ' CT< as time proceeds. 
Consequently, an initially localized probability distribution p(r, 0) oc 5 d (f) disperses as 

p{ft) = — / d d k exp{ik- f- \k\ a t) . (0 < a < 2). (9) 

{2K) a J 

This solution of the equation of motion fl7|) is called a Levy-stable distribution. As shown 
in Fig. [H the tails of such a Levy-stable distribution are more pronounced than those 
of an ordinary Gaussian distribution. Moreover, the width of the distribution (FS]) grows 
with time as t 1//<J , hence for a < 2 the dynamics of the random walk is superdiffusive. 

The aforementioned short-range cutoff in form of a minimal flight distance generates 
to lowest order an additional Laplacian in the equation of motion. For a < 2 this 
contribution k 2 is less relevant than |fc| CT and can be neglected asymptotically. For 
a > 2, however, the k 2 term becomes dominant and one obtains a Gaussian distribution 

= 77^m ex pf"^) • ( 10 ) 



(47Ti) rf / 2 *\ At J 

Thus, for a > o* = 2 anomalous diffusion crosses over to ordinary diffusion. 



2.2. Algebraically distributed waiting times 

In models for ordinary diffusion particles move spontaneously to nearest-neighbor sites at 
a constant rate, i.e., the time intervals At between subsequent moves are either constant 
or distributed exponentially like in Poissonian shot noise. In some situations, however, 
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Figure 1. Anomalous diffusion. Left: Trajectory of a random walk by Levy 
nights for a = 1.5. Right: Gaussian distribution (black) in comparison with 
Levy-stable distributions of the same width for a = 1.5 (red) and a = 1 (green). 



the distribution of waiting times in a random process is not exponential, instead one 
observes a broad distribution of waiting times. This may be described by a power-law 
distribution of the form 



where the exponent k is again a control parameter which plays a similar role as a in the 
case of spatial Levy nights. Interpreting the waiting times At as temporal displacements 
from event to event, they may be regarded as Levy flights in time. However, in contrast 
to the spatial Levy flights discussed above, temporal Levy flights have to be directed 
forward in time (At > 0) in order to respect causality. Such one-sided Levy flights differ 
from isotropic ones in that they are generated by a different type of fractional derivative 
which is non-symmetric under reflections. Its action on a Fourier mode is given by 



Note that in contrast to the spatial fractional derivative V CT this temporal operator brings 
down the nonsymmetric factor (iuj) K (rather than — \uj\ k ) in front of the exponential, 
reflecting the circumstance that temporal Levy flights are directed. For n > 1 the 
fractional derivative defined above crosses over to its short-range counterpart dt- The 
corresponding integral representation of the fractional derivative acting on a function 
git) reads 



P(At) ~ At 





(12) 



1 



/ 

Jo 



■oo 



dt'f 



I jJ-l-K 



[g(t) - g(t - 1')} 



(0 < K < 1) 



(13) 



where 




(14) 
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is the corresponding normalization constant. In this expression the function g(t) may 
be nonzero for all t 6 1, although in many applictions one uses the initial condition 
g{t) = for t < 0. 

The fractional derivative <9f can be used to write down an evolution equation for 
nearest-neighbor diffusion with waiting times 

d?p(f,t) = V 2 p(f,t). (15) 

It is important to note that such a random walk with algebraically distributed waiting 
times is no longer Markovian. This means that the temporal evolution is no longer 
determined by an initial configuration at a given instance of time, instead it depends 
on the whole history of the process. Therefore, as a convention troughout this paper, 
whenever we specify an initial condition at t = 0, we shall assume that the system was 
in the absorbing state for t < 0. 

For a localized initial condition at t = a formal solution of Eq. (|T5l) is given by 



p(f t) = - — — T / du / d d k—^ — '-. 16 

Remarkably this distribution is no longer normalized as time proceeds. In fact, the 
integral over the whole volume 

N(t)= fd d rp(r,t) = 1- rd, eXp( -^ 



2tt -{iuj 

r(i 



emr 1 (17) 

7T 

decays as N(t) ~ t K ~ l for k < 1. This is surprising since one expects N(t) to be 
conserved since the particle cannot disappear. However, for processes with waiting 
times the quantity p(f, t) no longer describes the overall density of all random walkers, 
instead it just gives the density of random walkers that are about to move at the next 
time step t, excluding the waiting ones. 

To explain this counterintuitive observation by an even simpler example let us 
consider the homogeneous equation 8* f(t) = 0. The function / may be thought of as a 
density of non-interacting particles transported forward in time by temporal Levy flights. 
Obviously, this equation is solved by f(t) = const, meaning that the density is conserved 
as time proceeds. For the corresponding Greens function, which is the solution of the 
inhomogeneous equation <9f f(t) = S(t), one would naively expect this conservation law 
to hold for t 7^ 0. However, the solution is given by f(uj) = (2ir)~ 1 / 2 (—iuj)~ K , i.e., 



/(*) 



const t K 1 for t > 
for t < 



which decays in time. This decay is a consequence of the non-Markovian properties of 
the fractional derivative. 
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2.3. Combination of spatial Levy flights and waiting times 



Combining spatial Levy flights and algebraically distributed waiting times one obtains 
an anomalous diffusion process with spatio-temporal long range interactions. In the 
continuum limit such a process is described by the partial differential equation [17, 18] 



d«p(r,t) = V°p{r,t) 



(19) 



A random walk according to this equation starting at the origin at t = is described 
by the formal solution 



p(f, t) 



1 



(2 



7T 



irf+1 



du / d a k 



d exp(ik ■ r — iut) 



\k\« 



(20) 



It is easy to show that this solution obeys the scaling form 



P(r, t) 



t 



1—K—dK/cr 



t K/a 



with the scaling function 

g{z) =g(z) 



(2< 



d+i 



du I d d k &Mik • *~ iu) 



2 d n ~(d+3)/2 



du 



dk 



\k\ a - (iuj) K 
d6 sin d - 2 9 



(21) 



(22) 



d _ 2 exp(ikz cos(9) — iu) 



k a 



[ico 



This expression is valid in the regime a < 2 and k < 1, where both fractional derivatives 
are relevant, while for a > 2 (or k > 1) the term k c (or (iuj) K ) has to be replaced by the 
corresponding short-range counterpart k 2 (or iu). Therefore, the dynamical exponent z, 
that relates spatial and temporal scales, is given by 



cr / 'hi if a < 2 and k < 1 

a if cr < 2 and K > 1 

2/k if cr > 2 and k < 1 

2 if a > 2 and k > 1 



(23) 



leading to the phase diagram for anomalous diffusion shown in Fig. [2j 



2.4- Numerical simulation of Levy flights 

In the following we describe how systems with long-range interactions can be 
implemented numerically. A Levy-distributed distance vector can be created as follows. 
At first one generates a normalized vector equally distributed on the unit sphere. 
This can be done efficiently without calling time-consuming trigonometric functions by 
repeatedly generating a random vector in a cube until its length is smaller than 1 and 
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Figure 2. Phase diagram of anomalous diffusion. Depending on the relevancy 
of the Levy operators, four different phases have to be distinguished. The 
background colors indicate whether the motion of the random walker is sub- or 
super diffusive. The dotted line marks the boundary of ballistic motion z = 1. 



then normalizing this vector (see code fragment in Tabled]). Next one has to generate 
a real- valued flight distance, i.e., a floating point number r distributed as 



by setting r := z^ 1 ^, where z is a random number drawn from a flat distribution 
between and 1. Note that r > 1 because 1 is the minimal flight distance. Finally the 
vector is multiplied by flight distance r. In continuum models the distance vector can 
be used directly while for models on a lattice it has to be rounded in such a way that 
the target position snaps in one of the lattice sites. 

Similarly, it is possible to generate algebraically distributed waiting times by setting 
At := z~ l ' K and, if necessary, rounding the result to a value matching the time 
increments of the dynamics. 

While generating Levy flights is simple, the real numerical challenge lies in the 
handling of finite size effects. Especially in those parts of the phase diagram, where the 
dynamics is superdiffusive, finite-size effects can be extremely dominant. There are two 
ways to master these difficulties: 

• The preferred method is to eliminate finite-size effects altogether by dynamical 
memory allocation. This method requires that the initial condition is localized, e.g. 
in form of a single particle at the origin. Instead of keeping a sparsely occupied 
lattice in form of a static array in the memory of the computer, the actual state 
of the system is stored in a dynamically generated list of particle coordinates. 
Suitable container classes are provided by various software packages such as the 




(24) 
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vector GenerateLevyFlight (double sigma) 
{ 

double x,y,z,n,r; 

do { x=rnd() ; y=rnd() ; z=rnd() ; } while (x*x+y*y+z*z>l) ; 

n = sqrt (x*x+y*y+z*z) ; 

r = pow(rnd() ,-1/sigma) ; 

return vector( x*r/n, y*r/n, z*r/n ); 

} 

Table 1. Code fragment in C++ which generates three-dimensional 
isotropically distributed Levy flights. It presupposes the existence of a suitable 
class for 3d vectors as well as a function rnd() that draws a floating-point 
random number from a flat distribution between and 1. The first three lines 
generate a normalized and isotropically distributed vector (x/n,y/n,z/n) on 
the unit sphere. This vector is then multiplied by the Levy flight distance r. 

C++-Standard Template Library (STL) [41]. On modern computers with 64 bit 
architecture, the effective lattice size of 2 Md possible coordinates can be considered 
as virtually infinite. 

• If a finite lattice is needed as e.g. in simulations starting with a homogeneous initial 
state, the proper handling of flight distances exceeding the system's boundaries 
is crucial. Comparing several possibilities it turned out that flights exceeding 
the boundaries should neither be discarded nor truncated to some maximal flight 
distance. Instead the cleanest scaling behavior is obtained if one applies periodic 
boundary conditions in their truest sense, namely, by wrapping a long-distance 
flight as many times around the system as needed to reach the target site. In one 
spatial dimension this is particularly simple since the coordinate of the target site 
is just i ± r modulo system size. 

2.5. Analyzing numerical simulations of systems with Levy flights 

Monte Carlo simulations of absorbing phase transitions are usually investigated by a 
palette of well-established numerical methods. For example, the value of the dynamical 
exponent z can be estimated by measuring the mean square distance from the origin 
(r 2 (t)) which is known to grow as t 2 l z . However, for systems with Levy-type interactions 
this method is no longer applicable. To see this consider a single Levy flight over 
distance r for a < 2. Although the flight distance r := z~ l l a of an individual move 
is a finite number, its second moment (r 2 ), i.e., the arithmetic average of the squared 
flight distance over an ensemble of many independent moves, is a diverging quantity. 
The same applies to the mean square displacement from the origin (r 2 (t)) after several 
subsequent moves. 

As a way out, it is useful to monitor the logarithm of the squared spreading 
distance and to compute the average e^ r2 ^ which amounts to study the geometric 
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Figure 3. Geometric and arithmetic average of the mean square distance from 
the origin r 2 (i) of an anomalously diffusing particle in 1+1 dimensions with 
a = 1.5. The left panel shows the numerically computed averages after 1000 
updates depending on the ensemble size M. As can be seen, the geometric 
average converges towards a well-defined limit whereas the arithmetic average 
does not. The right panel shows the average for various ensemble sizes as a 
function of time, demonstrating that the geometric average quickly reproduces 
the correct power law t 2//<T while the arithmetic average does not converge at 
all. 



instead of the arithmetic average. This average is finite and should grow as t 2 / z for the 
following reason. In systems exhibiting simple scaling (as opposed to multiscaling), the 
distribution of the mean square distance from the origin r(t) is expected to obey the 
scaling form 

P (r,t) = ^F(^) (25) 

with some scaling function F. This scaling form implies that the geometric average 

e (HrHt))) = exp jy d d r p^ j] ln r 2 1 ( 26 ) 

= exp jy & d rF{r) (21nr + -lnt)| oc exp[-lnt] 

(27) 

grows with time as t 2 ^ z which allows one to determine the exponent z. Therefore, we 
can use the geometric average in the very same way as the arithmetic one is used in the 
case of short-range models. 

Fig. [3] demonstrates the strikingly different behavior of arithmetic and geometric 
averages in the case of anomalous diffusion for a = 1.5. As can be seen, the arithmetic 
average of the squared distance from the origin does not converge as the size of the 
statistical ensemble increases whereas the geometric average converges quickly, showing 
the predicted scaling behavior. 
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Figure 4. Pair annihilation A + A — ► by Levy nights. The graph shows direct 
estimates and extrapolations for the decay exponent a as a function of a in 
comparison to the exact result (solid line). Data taken from Ref. [25]. 

3. Absorbing phase transitions with Levy-distributed interactions 

Having discussed anomalous diffusion of particles transported by Levy flights in 
combination with algebraically distributed waiting times, we now turn to systems where 
particles react with each other. The simplest one to start with is the pair annihilation 
process A + A — > (or, equivalently, the coagulation process A + A — > A). Later in this 
section we will consider more complicated systems with phase transitions into absorbing 
states. 

3.1. Pair annihilation by Levy flights 

In the ordinary annihilation process A + A — > with short-range interactions, the 
average particle density is known to decay as [42] 

( t- d l 2 for d < 2 , 
p{t) ~ \ t- 1 lnt for d = d c = 2 , (28) 
[ r 1 for d > 2 . 

This leads to the question how long-interactions by Levy flights change the dynamics, 
in particular the decay of the particle density and the value of the upper critical 
dimension. This problem was originally studied in Refs. [25,43], using both simulations 
and approximate theoretical techniques. 

Numerical simulation: 

The annihilation process A + A — > in 1+1 dimensions with Levy-distributed 
interactions may be realized on a linear chain of L sites with periodic boundary 
conditions, where each site i is either vacant (s{ = 0) or occupied by a particle (sj = 1). 
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The model evolves by random-sequential updates according to the following dynamical 
rules: 

(i) Select a random site i. If this site is vacant (sj = 0) discard the update, else proceed. 
Alternatively, the site may be randomly selected from a dynamically generated list 
of occupied sites which is more efficient. 

(ii) Draw a random number z e [0, 1] and generate a real-valued interaction distance 
r := z^ 1 ^. Determine the target site j := (i ± \r\) mod L in both directions 
with equal probability, where \r\ denotes truncation to the largest integer number 
smaller than r. 

(iii) If the target site j is empty, the particle diffuses from site i to site j, otherwise the 
two particles annihilate at site j. Both cases can be processed by setting Sj := 
and Sj := 1 — Sj. 

Performing numerical simulations one observes that the density decays algebraically as 

p(t) ~ r a((7) (29) 

with an exponent a that varies continuously with the control exponent a. Fig. H] shows 
numerical estimates of a(a) obtained in a simulation up to time t = 10 4 on a finite 
lattice with 2 16 sites averaged over at least 10 3 independent runs. It turns out that the 
expected power law (1291) is masked by considerable corrections to scaling, which requires 
to extrapolate the exponent to t — > oo. 



Field-theoretic analysis: 

On the field-theoretic level pair annihilation with Levy flights may be described by 
inserting the additional operator V 7 into the well-known field-theoretic action [42] 

S$, V] = J d d x dt {>ijj(dt - D N V 2 - D A V a )i) 

+ 2A# 2 + - poi>5(t) } , (30) 

where po is the initial homogeneous density at t — 0. Here the field is not simply 
related to the coarse-grained density field [44], although the average values of both 
fields coincide. Note that the corresponding action for the so-called coagulation process 
A + A — > A with Levy-distributed interactions would differ only in the coefficients 
of the reaction terms. Hence both processes are expected to have the same universal 
properties. 

For a < 2, power counting reveals that the upper critical dimension of the model 
is now reduced to d c = a < 2, hence for d > d c mean-field theory is expected to 
be quantitatively accurate, predicting the density decay p(t) ~ t~ x . Below the upper 
critical dimension, however, the renormalized reaction rate flows to an order e = o — d 
fixed point. Following the arguments of Ref. [42], for d < d c the only dimensionful 
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quantity is the time t, which, for o < 2, scales as [t] ~ r° ' . Hence, for a < 2, the density 



Note that for a > 2 the results cross over smoothly to the standard annihilation 
exponents ( l28i) . 

The field-theoretic prediction is shown in Fig. H] as a green solid line, which is in 
good agreement with the extrapolated numerical results. 

3.2. Directed Percolation with Levy flights in space 

Let us now turn to Directed Percolation (DP), which is the most important universality 
class of phase transitions into an absorbing state. DP is often used as a toy model 
for epidemic spreading: Infected individuals (active sites A) either infect one of their 
nearest neighbors (A — > 2 A) or they recover spontaneously (A — > 0). This competition 
between infection and recovery leads to a continuous phase transition out of equilibrium 
with universal properties. 

In more realistic situations, however, short-range interactions may be inadequate 
to describe the transport mechanism by which the infection is spread. For example, 
when an infectious disease is transported by insects, their motion is typically not a 
random walk, rather one observes occasional flights over long distances before the next 
infection occurs. Similar phenomena are expected when the spreading agent is carried 
by a turbulent flow. This motivated Mollison [12] to introduce a model for epidemic 
spreading by Levy flights, also called anomalous directed percolation. It turns out that 
Levy-distributed interactions do not destroy the transition, instead they may change its 
universal properties. 

Mean field approximation and phenomenological properties: 

Ordinary DP with short-range interactions is described by the Langevin equation [9] 



where V 2 is the Laplacian, p(f, t) is a coarse-grained density of infected individuals, and 
r, a, b, D are certain coefficients (r is introduced for later convenience and may be set 
to 1). The Gaussian noise £(r, t) accounts for density fluctuations. In order to ensure 
the existence of an absorbing state, the noise amplitude depends on the density of active 
sites. According to the central limit theorem the noise is described by the correlation 
function 



decays as 




(31) 



rd t p(f, t) = ap(f, t) - bp 2 (f, t) + DV 2 p(f, t) + f (r, t) 



(32) 



(£(f,t)Z(f,tf)) = 2cp(r,t) 5 d (f- ?)5(t - t') . 



(33) 
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Figure 5. Directed percolation with spatial Levy nights in 1+1 dimensions for 
various values of a. 



Following the strategy described above, long-range infections by Levy flights may be 
introduced by adding a term with a fractional derivative V' 7 (see Eq. 03)), leading to 
the modified Langevin equation 

rd tP = ap- bp 2 + DV 2 p + DV a p + £(f, t) (34) 



with the same noise as in Eq. (1321) . Dimensional analysis of Eq. ( 1341) immediately yields 
the mean field exponents 

oMF 1 MF -1 ,,MF 1 fnr\ 

p = 1 , u ± = a , z/|| = 1 , (35) 

and the upper critical dimension 

d c = 2a. (36) 

As in the previous case of anomalous pair annihilation, the mean field exponents are 
found to vary continuously with a and n. As expected, for a — > 2 and k — > 1 they cross 
over smoothly to the well-known mean field exponents of DP f3 = v\\= 2i>± = 1 in d > 4 
spatial dimensions. 

Spatial Levy-distributed interactions in a DP process have two different 
consequences. Firstly, the clusters are no longer connected, instead the long-range 
interactions allow any site to become active. This is demonstrated in Fig. [51 where 
the spatio-temporal evolution of a critical 1+1-dimensional process is shown for various 
values of a. Secondly, long-range interactions tend to enhance diffusive mixing, leading 
to a more homogeneous impression of the particle density as a decreases. Accordingly 
the upper critical dimension decreases since long-range interactions tend to diminish the 
influence of fluctuations. 

In the rightmost panel of Fig. [5] the process suddenly enters the absorbing state. 
Such finite-size effects are typical for phase transitions into absorbing states with long- 
range interactions. Within mean field theory, finite-size effects are expected to set at a 
typical time that scales as L a , where L is the lateral system size. Thus, especially for 
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small values of a finite-size effects can be extremely dominant. 



Field-theoretic renormalization group calculation: 

Below the upper critical dimension d < d c directed percolation with Levy flights is 
expected to be dominated by fluctuation effects. Here the critical exponents are no 
longer given by the mean field values (J35l) . instead they have to be computed by means 
of a field-theoretic renormalization group calculation. 

Starting point of a field-theoretic renormalization group calculation is an effective 
action S which is defined as the partition sum over all realizations of the density field 
if)(f,t) := p(r,t) and the noise field £(r, t), constrained to solutions of the Langevin 
equation (134|) . Integrating out the noise by introducing a Martin-Siggia- Rosen response 
field if)(r,t) and symmetrizing the coefficients of the cubic terms by rescaling the fields 
if) and if) one arrives at the field-theoretic action S[if>, ip] = S [ip, $\ + Si nt [ip, if)] consisting 
of a free contribution 

S [ip, $] = J d d x dt 4>{Td t -a- DV 2 - DV a )ip (37) 
and an interaction part 

S in thM = 9 f d d xdt $V 2 -$V) (38) 



with the coupling constant g = ybc. Note that the additional fractional derivative 
V CT appears only in the free part of the action. Therefore, the structure of the 
Feynman diagrams in a loop expansion is the same as in ordinary DP (see Fig. [6]), 
the only difference being that the free propagator in momentum space G^ p (k,u) = 
(Dk 2 — a — iruj)^ 1 is replaced by the generalized counterpart 

G (k,uj)= Dk 2 + Dk a - a - iruo . (39) 

Apart from this modification the loop integrals have exactly the same structure and the 
same multiplicities as in DP. In particular, the so-called rapidity reversal symmetry 

i>(k,oj) -> -uj) ${k,w) -> -ip(k, -^) (40) 

is still valid. As a consequence the fields if) and if) scale identically and the cubic vertices 
shown in Fig. [Hb-Eb renormalize exactly in the same way. 

As first shown by Janssen et al [24] , fractional derivatives used in interacting field 
theories have unexpected properties. Firstly, the loop corrections do not renormalize the 
Levy operator V CT , instead they generate short-range contributions, the most relevant 
being proportional to k 2 . This is the reason why the Laplacian has to be retained 
in the field-theoretic action for d < d c . In fact, even if this term was not included, 
it would be generated under renormalization group (RG) transformations. As shown 
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Figure 6. One-loop Feynman diagrams of a) the propagator and b)-c) the cubic 
vertices of directed percolation. The thin lines represent the free propagator 
Go(k,u) directed in time according to the arrow. The small circles denote the 
symmetrized bare cubic vertices ±g. Levy-distributed interactions modify the 
free propagator but they do not change the structure of the loop expansion. 



in the Appendix, the circumstance that the fractional derivative is not renormalized 
induces an additional scaling relation 



among the three standard exponents /3, u±, and v\\. The second observation by Janssen 
et al [24] is even more important. As in the previous examples of anomalous diffusion 
and annihilation, one would expect Levy-DP to cross over to ordinary DP at a = 2 
since this is the threshold from where on V CT behaves in the same way as a Laplacian 
V 2 . However, Janssen et al showed that this naive expectation is no longer correct in 
an interacting field theory: 

"Because the anomalous susceptibility exponent, the analog of the Fisher exponent in 
critical equilibrium phenomena, is negative here, the continuous crossover between 
long-range and short-range behavior arises at a Levy exponent [a] greater than 2. " 

This means that in an interacting field theory the Levy operator V CT does not cross over 
to V 2 at a = 2, rather its long-range nature prevails even beyond this threshold up to 
some a* > 2. Assuming a continuous crossover of the critical exponents the threshold a* 
can be computed by plugging the numerically known critical exponents of short-range 
DP into the scaling relation (I4ip and to solve it for a, giving 



with a* = 2.0766(2) in one, cr* ~ 2.2 in two, and a* = 2 + e/12 in d = 4 — e spatial 
dimensions. 

Numerical results: 

As shown in [24] (see also Appendix) a field-theoretical renormalization group 




(41) 



a* = d + z - 2f3/v ± 



(42) 
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Figure 7. Numerically determined critical exponents of directed percolation 
with spatial Levy flights. The solid lines indicate the mean field predictions 
and the field-theoretic results to one- loop order (data taken from [25]). 



calculation to one-loop order for DP with Levy flights yields the critical exponents 

= 1-^ + (e 2 ), v±= l+*L +0 (?), ,,, = 1 + ^ + (e 2 ) .(43) 

In order to verify these predictions numerically, it is advisable to perform epidemic 
simulations starting with a single seed [45] using a dynamically generated list of active 
sites. As described in Sect. 12.41 this simulation technique virtually eliminates finite size 
effects, which would be very strong especially when a is small. As usual in this type 
of simulations, one measures the survival probability P s (t) ~ t~ s , the average number 
of particles N(t) ~ t e , and the mean square spreading from the origin R 2 (t) ~ t 2 / z . 
Since the latter diverges in the presence of Levy flights, the arithmetic average has to 
be replaced by a geometric average, as explained in Sect. 12.51 Having determined the 
critical point for each value of a separately, the power law decay of these quantities 
allows one to determine the exponents a — S — z — u»/u±, and 6 = - — 2(3/ v\\. 

In Fig. [7] the analytical predictions are compared with the numerical estimates of 
the critical exponents in 1+1 dimensions. For o < 1/2 very long Levy flights are so 
dominant that they wash out any fluctuations, hence the system is effectively described 
by mean field exponents. In an intermediate regime 1/2 < a < a* the exponents 
vary continuously. Here the field-theoretic one-loop results correspond to the tangents 
(solid lines) at the left boundary, being in qualitative agreement with the numerical 
findings. Finally, for o > a* the numerical results confirm the crossover to ordinary DP, 
although the results are not precise enough to verify the predicted value of a*. This 
example shows that the strength of field theory is not in a quantitative prediction of 
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the exponents, rather the field-theoretic approach provides a conceptual understanding 
of the underlying process, allows one to derive exact scaling relations, and explains the 
origin of universality. 



3.3. Directed Percolation with waiting times 



As in the case of pure diffusion, it is straight-forward to generalize ordinary DP as to 
include algebraically distributed waiting times. The corresponding Langevin equation 
has now a fractional derivative in time: 



8?p + d tP = ap- bp 2 + DV 2 p + £(f, t) 



(44) 



As before, we keep the operator d t p which reflects the short-range cutoff of the 
distribution. In analogy to spatial Levy flights, it is also natural to expect that 
such a term will be generated under renormalization group transformations. Ignoring 
fluctuation effects, dimensional analysis yields the mean field exponents 



P 



MF 



V 



MF 
_L 



1/2. 



MF 

'II 



(45) 



and the upper critical dimension d c = 6—2/k. Note that it is now the temporal exponent 
v\\ that varies continuously with k. Like spatial Levy flights, algebraically distributed 
waiting times decrease the upper critical dimension, meaning that they tend to wash 
out fluctuation effects. 

Using similar arguments as in the case of spatial Levy-distributed interactions one 
can show that the fractional operator d£ does not renormalize itself, instead it renders 
corrections of the form iuo and hence renormalizes the short-range operator dt- This 
leads to the scaling relation analogous to Eq. (14TI) : 



Un(l-K)+dv ± - 2/3 = 0. 



(46) 



Recently Jimenez-Dalmaroni [30] presented a one-loop renormalzation group calculation, 
reporting that the critical exponents in d = d c — e dimensions are given by 







1 + 
1 

2 + 

1 

K 



\{2k 



tan(7w/2) 
8(2k - 1) 
Ktan(^K) 

8(2«- r 



where 



F 1 (k) 



k 2 Fz(k) tan(f/t) 



, )(2 1\ sin(TTK) 
K/sin 2 (|/t) 



[F x {k) + F 2 («)] 
Fi(«) -F 2 {k 
0{e 2 ). 



3k 



2 

e 



2F, k 



kF 3 (k) 
4 

+ 0(e 2 ). 



kF z (k i 



+ 0(e 2 



(47) 



/7T 
COS ^ — AC 



Ak 2 - 3k + 1 



3k- 1 



COS 



3tt 



-K 



(4£ 



F 3 (k) = (32 - d c )F x {K) + d c F 2 (K). 
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3.4- Directed Percolation with spatio-temporal long-range interactions 



Finally it is possible to consider the mixed case of epidemic spreading by spatial Levy 
flights combined with algebraically distributed waiting times. Power counting in the 
corresponding Langevin equation 

fd?p = ap- bp 2 + DV a p + f (r, t) (49) 

with the same noise as in Eq. ( |33i) yields the mean field exponents 

oMF -i MF -1 MF , -1 /r n \ 

p = 1 , u ± = a , u» — K (5U) 
and the upper critical dimension 

d c = 3a - - . (51) 

In particular, the dynamical exponent locks onto the ratio z = o~/k. 

A field-theoretic approach to the mixed case can be found in [31]. It turns out that 
in this case both scaling relations (]4"Ij) and (|46p are simultaneously valid, hence only one 
of the three exponents is independent. To one- loop order one finds the results 

(3=l- e - + 0(e 2 ), ^|| = l + ^ + 0(e 2 ), u ± = I + ± + 0(e 2 ). (52) 

The phase diagram of a 1+1-dimensional DP process with spatio-temporal long-range 
interactions controlled by o and k is sketched in Fig. [HJ It comprises seven different 
phases: 

(i) a mixed mean field phase for small a and k which is governed both by Levy flights 
and waiting times. In this phase fluctuation effects are irrelevant and the exponents 
are given by Eq. (|50|) . 

(ii) a mean field phase for k > 1 governed solely by spatial Levy flights. In this regime 
waiting times described by the operator d* are irrelevant compared to d t and the 
exponents are given by Eq. ( |35l) . 

(iii) a mean field phase for a > 2 governed solely by waiting times. Here spatial Levy 
flights associated with the operator V CT become irrelevant and the exponents are 
given by Eq. (j4"5"l) . 

(iv) a non-trivial fluctuating mixed phase in the center, where spatial Levy flights and 
incubation times are both relevant. Here both scaling relations ( HIT) and ( 1461) are 
simultaneously valid, the upper critical dimension is d c = 3a — -, and the exponents 
are given to one-loop order by Eq. (l52~l) . 

(v) a non-trivial fluctuating phase dominated exclusively by spatial Levy flights, where 
incubation times are still irrelevant so that the results of Refs. [24, 25] apply, see 
Eqs. flU}, flUD and fl43]). 
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Figure 8. Phase diagram for epidemic spreading with spatio-temporal long- 
range interactions in one dimension. The seven phases and their boundaries 
are explained in the text. 



(vi) a non-trivial fluctuating phase dominated by incubation times, where spatial Levy 
nights are irrelevant so that the results of Ref. [30] apply, see above. 

(vii) a short-range regime where the process behaves effectively in the same way as an 
ordinary DP process. In 1+1 dimensions the corresponding exponents are given by 
f3 DP ~ 0.2765, z/^ p ~ 1.097, and z/ DP ~ 1.734. 



3.5. Synchronization of extended chaotic systems with Levy-type interactions 

Long-range interactions by Levy nights can also be studied in various other contexts. 
For example, in a recent paper Tessone et al [46] studied a coupled map lattice defined 
by the update rule 

Xi (t + 1) = (1 - 7 ) F(xi(t)) + "fFifat)) (53) 
Vi {t + 1) = (1 - 7 ) F(Ut)) + lF( Xi (t)) . 

Here t and i are discrete temporal and spatial indices while Xi(t),yi(t) are the state 
variables of two linear chains of coupled maps with periodic boundary conditions. The 
map F is chaotic; a typical choice would be the Bernoulli map F(x) = (2a;) modi. In 
addition, the two chains are coupled crosswise controlled by the parameter 7. Depending 
on 7 the two chains either synchronize or they evolve into independent chaotic states. 
The transition from independent dynamics to synchronization is characterized by the 
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Hamming distance and can be interpreted as a transition from a fluctuating phase into 
an absorbing state. 

In the update rules the functions Xi(t) and yi(t) are linear combinations of the state 
variables and are used to encode different coupling schemes. For example, ordinary 
short-range interactions would correspond to Zi(t) = (1 — e)zj + e(zj_i + z i+1 ). Obviously 
the generalized coupling scheme for non-local Levy-distributed interactions would be 
given by 

m = (i-^ + JLjEii+jtii, (54) 

where j is the interaction distance with a suitable cutoff at the system size and i](a) is the 
corresponding normalization. However, using the language of DP, this update algorithm 
first selects a non-infected site and looks from where an infection could come, thus 
summing over all sites of the chain. Consequently, each of these updates requires 0(L 2 ) 
operations which strongly limits the numerically manageable system sizes. Contrarily, 
the update rules discussed in the previous section determine the target site to where the 
infection goes which takes only O(L) operations. 

To circumvent this problem Tessone et al present a remarkable approach which 
could be pathbreaking for many similar problems. The idea is to replace the fully- 
coupled interaction scheme of Eq. (1541) by a diluted coupling scheme which has 
nethertheless the same scaling properties. More specifically, Tessone et al restrict the 
couplings to a sparse set of allowed interaction ranges 

3m(q) = q m -l, (55) 

where q is a parameter with values typically chosen as q = 2,4, and 8. The index m 
enumerates the allowed interaction distances and runs from 1 to M = log 9 (L/2). The 
coupling scheme is then given by 

Al 

which differs from the fully-coupled one in Eq. (|54j) in so far as the interaction distance 
in the denominator is raised to the power a instead of 1 + a. As shown Tessone et 
al. this modification ensures that the diluted coupling scheme has effectively the same 
statistical weight as in the fully-coupled case. Because of the dilution the computational 
effort per site reduces to 0{L log L) operations per update which allows one to simulate 
much larger systems. 

Tessone et al studied the scaling behavior of the synchronization error Wi(t) = 
\xi(t) — Ui{t)\ in simulations of a coupled map lattice with up to 4 x 10 6 maps (sites). 
Comparing their numerical results with those of Ref. [25] they find that the system 
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belongs to the same universality class as anomalous directed percolation with spatial 
Levy-distributed interactions. 

3.6. Parity- conserving processes with spatial Levy flights 

Although so far most papers focused on generalizations of directed percolation, the 
concept of long-range interactions can be applied to any system with a phase transition 
into absorbing states. Another important universality class of such phase transitions 
is the so-called parity-conserving (PC) class [47-50] in 1+1 dimensions, sometimes also 
referred to as the voter universality class [51]. This process is described by the reaction- 
diffusion scheme 

2A -> 3 A , 2A -> (57) 

and exhibits a non-equilibrium phase transition from a fluctuating phase with a 
stationary density of particles to an absorbing phase dominated by pair annihilation with 
an algebraic decay ~ t~ x l 2 . The upper critical dimension of this process, above which 
mean-field behavior prevails, is d c — 2, while in an intermediate regime d' c < d < d c 
with d' c — 4/3 the dynamics is governed by the annihilation fixed point [42] with a zero 
branching rate at the critical point, see Eq. (1281) . Only below d' c the critical branching 
rate becomes nonzero and the transition is characterized by non-trivial properties. 
Unfortunately, this case cannot be described by a perturbative renormalization group 
calculation [49]. 

It is now straight-forward to replace nearest-neighbor couplings by non-local 
interactions according to a symmetric Levy distribution controlled by a. In a field- 
theoretical and numerical study in 1 + 1 dimensions Vernon and Howard [52] observed 
the following phase structure (see Fig. [9j): 

(a) For a > 2.5 the process is effectively a short-ranged and the transition belongs to 
the usual PC class. 

(b) In a range 1.5 < a < 2.5 one observes a nontrivial phase with continuously varying 
exponents. 

(c) For a < 1.5 the critical branching rate is zero and one finds the same behavior as 
for pair annihilation with Levy flights as described in Sect. 13.11 where the density 
decay according to Eq. (l3"Tj) . 

(d) For a < 1 the process displays mean field behavior with the density decaying away 
as p ~ 1/t at zero branching rate. 

The observed phenomenology in 1+1 dimensions suggests the general phase 
structure shown in Fig. [9] Vernon and Howard show that the phase boundary, where 
the branching rate becomes zero, is given by cr' c (d) = 3d/2, while the precise form of the 
phase boundary between short- and long-ranged branching annihilating random walks 
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Figure 9. General phase structure of parity-conserving branching annihilating 
random walks with Levy-distributed long range interactions. 



(the dashed line in Fig. [9]) is not yet known. However, taking field-theoretic action given 
in Ref. [52] and assuming that the fractional derivative is not renormalized, it is near at 
hand to conjecture the scaling form 

a = d + z(l-a-6) (58) 

with 5 = 0. This relation should hold in the phases (b) and (c), where the exponents 
vary continuously with a, and gives the correct phase boundary of the annihilation- 
dominated regime. If this relation were correct, it would give the phase boundary of 
phase (b) (the dashed line) by plugging in the numerically known exponents of the PC 
class and solving the equation for a. In one dimension this would give the threshold 
er c m 2.50, which is in qualitative agreement with the numerical simulations of Ref. [52]. 



4. Restricted long-range interactions 

Recently phase transitions into absorbing states with long-range interactions emerged 
in a different context, namely, in wetting processes out of equilibrium [33-39]. Here 
one studies an interface growing on a ^-dimensional hard-core substrate. The interface 
evolves in time by deposition and evaporation processes which generally violate detailed 
balance. Depending on their rates the interface is either pinned to the substrate or it 
detaches, leading to a so-called wetting transition. 

A schematic typical configuration of such an interface close to the wetting transition 
is shown in Fig. [TQ1 It consists of detached domains separated by pinned segments. The 
dynamics of the interface is such that each detached domain can shrink or expand from 
its edges. New domains may be created by an unbinding process of bound sites, and 
two or more detached domains may merge into a single larger one. 
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Figure 10. Non-equilibrium wetting process close to the unbinding transition. 
The contact points, where the interface touches the substrate (red points), may 
be interpreted as active sites of a directed percolation process. Figure taken 
from [37]. 



In principle a detached domain may also split into two domains when an unbound 
segment of the interface touches the substrate at the interior of the domain (marked by 
the arrow in Fig. [TUl) . Let us first consider those regions of the parameter space where 
such splitting processes are highly suppressed. For example, this is the case when an 
unbound interface tends to move away from the substrate while it is held bound to it 
by some short range attractive interaction. In this case the dynamics of the interface 
may very well be described by a directed percolation process in which the active sites 
are those bound to the substrate. In fact, in some models [53-55] it was observed that 
the depinning transition belongs to the DP universality class. However, the observation 
of first-order transitions under similar conditions [34] suggests the existence of effective 
long range interactions between the pinning sites, which may provide a mechanism 
that can lead to first-order transitions in one-dimensional models [56]. To explain this 
mechanism, Ginelli et al [37] introduced a variant of directed percolation with long- 
range interactions called a-process, which will be described below. 

On the other hand, in those regions of the parameter space where the splitting 
of unbound segments by reattachment in the middle is likely, a different long-range 
mechanism is needed to describe the dynamics of the contact points. On the level of a 
DP-like process such long-range interactions can be modeled by truncated Levy nights, 
called a-process, which will be discussed in Sect. 14.21 

Both generalizations, the a- and the a-process, differ significantly from DP with 
Levy-distributed interactions and waiting times discussed in the previous section. As 
will be shown, they give rise to fractional powers of the density field (instead of fractional 
derivatives) in the corresponding Langevin equation. 
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4-1. Local transport with non-local rate dependence: The a -process 

As ordinary DP, the a-process involves spontaneous recovery 1 — >• and infection of 
nearest neighbors 10/01 — > 11 and thus can be considered as a local process. However, in 
contrast to DP, where the rates for infection and recovery are constant, the infection rate 
now depends on the actual distance to the nearest neighbor, i.e., it varies with the size of 
the unbound segment in the corresponding unbinding transition, introducing an effective 
long-range interaction. More specifically, the infection process may be implemented as 

1 — ► with rate 1 , 

e l [lO 1 ] -> 0^11 [HO*- 1 ] with rate , 

where X(£) = A(l + a/£ a ) depends algebraically on the size I of an inactive island of 
size £. The process is controlled by three parameters, namely, the mean infection rate A, 
the control exponent a, and an amplitude a controlling the intensity of the long-range 
interactions. For a = the model reduces to the ordinary contact process with a tran- 
sition at A c = 1.64892(8) which belongs to the DP universality class [1]. 

Mean field approximation: 

Denoting by p(t) the average density of active sites at time t, the mean field equation 
in the thermodynamic limit L — > oo takes the form 

| = (A-l) P -Ap 2 + Aap 2 f;ii-^. (59) 

e=i 

For a > 1 the sum on the r.h.s. is finite in the limit p — > and its contribution amounts 
to a renormalization of the coefficient of the p 2 term so that the mean field behavior 
of a standard DP process with short range interaction is recovered, provided that a is 
small enough. On the other hand, for < a < 1, the leading contribution of the sum 
can be captured by replacing it with an integral over dl which, to leading order in p, 
reduces to 

= (A - l)p + aAr(l - a)p 1+a - Ap 2 . (60) 

Here T is the standard Gamma function. The leading nonlinear term in this equation 
involves a non-integer power of the density field, reflecting the long range nature of 
the interactions. Since its coefficient is positive, Eq. (1601 does not admit fixed point 
solutions for arbitrarily small densities, so that the transition to the absorbing state 
becomes discontinuous. 

Numerical simulation 

Numerical simulations suggest that the mean field prediction (DP for a > 1, first-order 
transition for a < 1) is qualitatively correct even in the full model in 1+1 dimensions [37]. 
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Figure 11. Snapshots of the <r-process close to criticality starting with a single 
seed for various values of a. 

To detect the crossover from a continuous to a discontinuous transition in one spatial 
dimension, it is useful to perform simulations starting with a single seed and to generate 
a cluster of active sites at criticality. As usual, DP behavior manifests itself by sparsely 
occupied clusters. Contrarily, first-order transitions can be recognized by the emergence 
of a single compact domain. As there is no surface tension in one- dimensional systems, 
these domains fluctuate without bias, hence the process becomes equivalent to a process 
called compact directed percolation [57]. In fact, as illustrated in Fig. |TTl for a < 1 the 
generated clusters become compact on large scales. 

Numerically it is easy to identify these two phases since the exponents of DP in 
1 + 1 dimensions (a ~ 0.159, 9 « 0.313, and z « 1.58) differ significantly from those of 
compact directed percolation (a = 1/2, 9 = 0, and z = 2). The simulations carried out 
in Ref. [37] clearly indicate that the crossover takes place at a close to 1. 

4-2. Truncated Levy flights: The a -process 

If the parameters of the wetting process are chosen in such a way that an unbound 
segment of the interface may reattach to the substrate in the interior (marked by the 
arrow in Fig.fTOl). an effective model for the dynamics of binding sites at the bottom layer 
has to incorporate an additional type of long-range interactions. In numerical studies of 
wetting transitions, the probability of reattachment at the critical point of the unbinding 
transition was found to decay algebraically with the distance from the nearest binding 
site. This motivated Ginelli et al to study a model called a-process which mimics these 
reattachments [58]. The model evolves by random sequential updates with the transition 
rates 

l->0 with rate 1, (61) 
— > 1 with rate q/l a . 
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The rate at which such a process takes place decays algebraically with the distance I 
from the nearest active site (measured in lattice site units), reflecting the long range 
nature of the interaction. Note that in contrast to the cr-model, where infections are 
local (although taking place with a non-locally determined rate), the a process admits 
infections over long distances. Obviously, it reduces to ordinary DP in the limit a — > oo. 

It should be emphasized that the a-process differs significantly from spreading pro- 
cesses with Levy-distributed interactions discussed in the previous section. Since an 
inactive site can only be activated by the nearest active site, the interaction range is 
effectively cut off by half the actual size of the corresponding island of inactive sites. 
This is to say that infections in the a-process are mediated by truncated Levy flights, 
which cannot overtake neighboring active sites. 

Mean field approximation: 

Assuming different sites to be uncorrelated, the probability of an inactive site to lie at 
a distance I from the nearest active site is (1 — (1 — p) 2 )(l — p) 2l ~ 2 - Summing up the 
contributions for all distances I, the mean-field equation governing the dynamics of the 
density of active sites takes the form 

d t p(t) = -pit) + qp(t)(2 - p(t)) ~ -J^— ■ ( 62 ) 

1=1 1 

Turning the sum into an integral the leading terms in p(t) in this equation can be 
calculated. It is found that for a < 1 and, to leading order in p(t), the mean field 
equation turns into d t p(t) = rp[t) + up(t) a , where r and u are constants with u > 0. 
Since the leading term in this equation is p a , and since its coefficient is positive, the 
absorbing state (p = 0) is always unstable and no transition takes place for any finite 
q. For 1 < a < 2, however, the leading term in Eq. (1621) is the linear one and the 
coefficient u is negative. This results in a continuous transition different from DP which 
takes place at r = corresponding to a non- vanishing rate q. 

The above mean field equation has to be extended by a term that describes the 
transport mechanism. In systems with non-conserved short-range interactions this term 
is a Laplacian, while in systems with unrestricted Levy-distributed interactions it would 
be a fractional derivative (see Eq. (jSJ)). In the a-process, however, the interactions are 
power-law distributed but their range is cut off at the distance to the closest active 
site. Therefore, in the large-scale limit (which the mean field equations are supposed to 
describe) the interactions are effectively short-ranged and thus described by a Laplacian. 
Instead the effect of the non-local infection process is to make the diffusion coefficient 
D dependent on the density p. As shown in Ref. [58] , for 1 < a < 2 this reasoning leads 
to the mean field equation 

d t p{f, t) = rp(r, t) + u[p{f, t)} a + [p(f, t)} a - 3 V 2 p(f, t) (63) 
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for which dimensional analysis yields the complete set of mean-field exponents 

< F = [3 MF = — *— , vr = 1 , z MF = a - 1 . (64) 
a — 1 11 

On the other hand, for 2 < a < 3 the leading non-linear term is the quadratic one, 
leading to the lowest-order mean field equation 

d tP (f, t) = rp(r, t) + v[p(f, t)] 2 + [p(f, t)] a - 3 V 2 p(r, t). (65) 
with the critical exponents 

z/ mf = 2 _« /3 mp = i/ mf = 1> 2 MF = 2/(4 - a) . (66) 

Finally, for a > 3 the leading-order mean field equation is the same as the one of DP, 
hence we expect the process to cross over to ordinary directed percolation with z — 2. 
This phase structure is also observed numerically, although with different values of the 
critical exponents [58]. 



Combination of the a- and the a process: 

Both models, the a-process and the a-process, describe different aspects of the effective 
dynamics of binding sites in a non-equilibrium wetting process. In most situations, 
however, one expects a mixture of both mechanisms. Such a combined a — a process 
may be modeled by the transitions 

1 ^0 with rate 1 

10,01 11 with rate q(l + a/rrf ) (67) 

000 -> 010 with rate bq/l a 

where m is the size of the inactive island and I denotes the distance from the nearest 
active site. Here q is the control parameter, a and a are the control exponents, and a 
and b are constants. The "pure" a- and cr-processes are recovered in the limits a — > oo 
and a — > oo, respectively. 

Within a mean field approach, this combined a-a process is described by the 
equation 

00 00 (~\ \2l— 1 

dtp = -P + gp 2 £(1 + ^)(1 - P) m + bq P {2 -p)J2 [ ~j P J -(68) 

lib V 

ra=l 1=2 

In the regime of interest, namely for < a < 1 and 1 < a < 2, the mean-field equation 
can be written to leading order as 

d t p = rp + pp 1+a - up a + 0(p 2 ) , (69) 



where 



22— a 2° ^ 

q-l + bq -, p = aqT(l-<T), u = bq -V(2-a) (70) 

a — 1 a — 1 
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Figure 12. Phase diagram of the combined a-a process 



are g-dependent constants. As shown in Ref . [58] , the corresponding mean field diagram 
(see Fig. [121) contains four phases, namely, 

(i) a DP phase for a > 1, a > 3, 

(ii) a phase for o > a — 1 and 1 < a < 2, where the transition is second order with 
continuously varying exponents. 

(iii) a phase for o > 1 and 2 < a < 3, where the transition is second order but where z 
is the only continuously varying exponent. 

(iv) a phase for o < a — 1 and < a < 1, where the transition is discontinuous. 
5. Concluding remarks 

Continuous phase transitions into absorbing states continue to fascinate theoretical 
physicists because of their universal properties which are determined by symmetries 
and conservation laws rather than microscopic details. On a field-theoretic level the 
origin of universality can be understood in terms of relevant and irrelevant operators in 
the corresponding action. Usually it is assumed that the process evolves by a Markovian 
dynamics with local interactions, as described to leading order by differential operators 
such as d t and V 2 . Generalizing these models by introducing long-range interactions 
and non-Markovian effects with power-law characteristics these universal properties may 
change. 

In Sect. [3] we have discussed how unrestricted long-range interactions over a 
distance r distributed as P(r) ~ r ~d-<r an j wa iting times At with a power-law 
distribution -P(At) ~ affect phase transition belonging to the universality class 

of directed percolation. Generically one observes the following phenomenology: If the 
control exponent o (or k) is larger than a certain upper threshold cr* (or k*) the model 
behaves effectively as if the interactions where short-ranged. If the exponent is smaller 
than a certain lower threshold a c (or k c ), the couplings become so long-ranged that 
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an effective mean-field behavior is observed. In between there is a non-trivial phase in 
which the critical exponents (3, u±, and uu vary continuously with the control exponent. 
Generally long-range interactions in space as well as in time tend to reduce fluctuation 
effects, thereby decreasing the upper critical dimension of the process. 

On a field-theoretical level Levy-distributed long range interactions can be studied 
by adding terms with fractional derivatives. Generalizing directed percolation in this 
way one obtains a renormalizable field theory describing a spectrum of universality 
classes parameterized by a and k. Moreover, it is found that the additional terms 
involving fractional derivatives are not renormalized, implying exact scaling relations 
among the critical exponents in the non-trivial phase. 

These scaling relations can be used to compute the upper thresholds a* and k*, 
where the crossover to the short-range behavior takes place. Surprisingly one finds that 
a* > 2 and k* > 1. This example demonstrates that in an interacting field theory the 
thresholds of a and k, where the fractional derivatives V 7 and 8* cross over to their 
short-range counterparts V 2 and d t , may be shifted away from the naive thresholds 
a = 2 and k = 1 that would be obtained if these operators were acting on smooth 
functions. More recently the same issue was debated in the context of the Ising model 
with Levy-distributed couplings [59-61] without being aware that the analogous problem 
was already solved for DP. 

In Sect. H] we have discussed a different type of long range interactions motivated by 
recent studies of non-equilibrium wetting processes. In the first case, called cr-process, 
the interactions are local but the coupling constant has a non-local dependence. This 
modification induces a crossover from directed percolation to a first-order transition. 
In the second case, called a-process, the interactions are non-local and power-law 
distributed but in contrast to the unrestricted interactions discussed in Sect. |3]they are 
now cut off at the actual distance to the nearest particle. It turns out that fractional 
derivatives are no longer suitable to describe such truncated Levy flights, instead one 
gets fractional powers of the density field and an algebraically varying diffusion constant 
in the corresponding Langevin equation. 

Obviously, the study of react ion- diffusion models and non-equilibrium critical 
phenomena with long-range interactions is still at its beginning. As demonstrated in 
this review, such problems are very fascinating as they can be treated both by numerical 
and analytical methods. Since long-range transport mechanisms are frequently observed 
in complex systems, it is the hope that the theoretical insight from simple models will 
be useful for the understanding of experiments in the future. 
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Appendix A. 



Field-theoretic approach to directed percolation with Levy flights 



This appendix sketches the field-theoretic approach to directed percolation with spatial 
Levy flights. This problem was first solved by Janssen et. al. [24] using dimensional 
regularization techniques. Here we present an alternative calculation using Wilsons 
momentum shell approach, leading to the same results. 

Starting point is the field-theoretic action 

5[^,V]= Jd d rdt^(Td t -a-DV 2 - DV a )ip (A.l) 

+ giH>{il> - i>) - po5{t)$ 



Here r is a coefficient fixing the time scale, a plays the role of the reduced percolation 
probability, D and D are the coefficients for short-range diffusion and Levy flights, 
g is the symmetrized nonlinear coupling constant, and po denotes the initial density 
of particles at t — 0. The only properties of the fractional derivative V CT used in the 
following renormalization group calculation are its scaling dimension as well as its action 
in momentum space 

VV^ = -|A;|V^\ (A.2) 

Scaling transformation: 

The first step of Wilsons renormalization group procedure is to rescale the action by 
f^bf, t^b z t, p ^^ PO Po, V->&~ d *V'> ^p^b'^^p, (A.3) 

where b is a dilatation parameter, z = v\\/v±_ is the dynamical critical exponent, while d^, 
d$, and d Po are scaling dimensions associated with the fields. Because of the well-known 
rapidity reversal symmetry of DP, which holds also in the presence of Levy flights, the 
fields ip and ip scale identically, meaning that d^ = d^. It is convenient to express these 
scaling dimensions in terms of the anomalous dimensions rj^ and r] Po by 

d^ = ^{d + v^), d po = ^(d + rj po ) . (A.4) 



The scaling transformation flA.31) leads to a change of the coefficients in the action flA.lj) 



For an infinitesimal dilatation b = 1 + i these changes are given to lowest order by 

d t hiT = 

di In a = z — rj^ 

dnhiD = z — rj^j — 2 

dg In D = z — rjip — a (A. 5) 
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Figure Al. Schematic sketch of Wilsons renormalization group procedure. 
A minimal length representing the lattice spacing is introduced by imposing 
an upper cutoff A in momentum space. In a first step all quantities are 
rescaled infinitesimally, modifying the cutoff scale. In order to compensate 
this modification, the high-wavenumber modes within the momentum shell are 
integrated out in a second step and the resulting contributions are absorbed to 
lowest order in the coefficients of the action. 



delng — z — 3^/2 — d/2 
dilnpo = - (Vip + Vpo)/ 2 

These are the renormalization group equations on tree level which are expected to be 
valid above the upper critical dimension d c = 2a, where loop corrections are irrelevant. 
The fixed point of these equations is z* = 2 and rf^ = r]* o = 0. Together with (3 MF = 1 
and = (3/u± they determine the mean field exponents 

oMF i ,,MF -1 ,,MF 1 /a (t\ 

(J = 1 , U ± —a , Uu =1 (A-Dj 

which are expected to be exact above the upper critical dimension d c = 2a. Moreover, 
it turns out that the critical initial slip exponent 

9 = 7 lp0 - m , (A.7) 

2z ' K ' 

which describes the average increase of the particle number in simulations starting with 
a single seed, vanishes at tree level. 



One-loop diagrams: 

Wilsons renormalization group approach computes loop corrections by introducing an 
upper cutoff \k\ < A in momentum space. Since the rescaling procedure changes the 
cutoff by A — > 6 _1 A a second step is required that sends the cutoff back to its original 
value (see Fig. |A1[) . This can be done integrating out all fast oscillating modes within 
the (infinitesimal) momentum shell A < \k\ < b~ 1 A and then to absorb the resulting 
contributions to lowest order into the coefficients of the action. 

In Wilsons renormalization group approach to ordinary DP (see e.g. Ref. [62]) the 
propagator in momentum space is renormalized to one-loop order by 



-G \k,uy en - = -G \k,u) - 




(A.8) 
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where 



G = (Dk 2 -a-iruj)- 1 (A.9) 
is the free propagator and I p denotes the loop integral 



^H/^^ + ^W)^-^ -.<). (A-10) 

Here the symbol '>' indicates that the integration is restricted to the momentum shell 
A < \k'\ < A/b. Similarly, the non-linear vertices are renormalized by 

-g Ten =-g + 2g 3 I v (A.ll) 
with the loop integral 

j-^Gl(k,u)G (-k,-u). (A.12) 

As discussed in the main text, spatial Levy flights can be introduced by replacing the 
free propagator (1A.9I) with 

G = (Dk 2 + D\k\ a - a - iruj)- 1 (A.13) 

while the structure of the loop integrals and their multiplicities do not change. There- 
fore, all what has to be done is to re-evaluate the loop integrals (lA.lOj) and (1A.12|) . As 
explained above, it is important to retain the short-range term k 2 in the action since 
the Levy operator does not renormalize itself, instead it renormalizes the Laplacian. 

Calculation of the loop integrals 

Let us first evaluate the propagator loop integral (jA.lOp . With the abbreviation 



A ± = D{k' ± k/2) 2 + D\k' ± k/2\ a - a - iru/2 (A.14) 

this integral can be expressed as 

„- . f d d k'du> 1 
Ip{k,u) 



(2TT) d+1 (A+ - ITU'){A- - ITU') 

1 f ddk> IT 1 a- - ( A - 15 ) 



(2n) d r y> A+ + A 
Since k' is much larger than k we may expand the denominator by 

i£'+~r+i£'-~r= (A.i6) 

= 2\kT + -\kv- 2 k 2 + -(- - 1) ik'i^k 2 cos 2 e + o(k A ) 

4 2 y 2 ! 1 1 y 1 
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where is the angle between k and k' . The loop integral over the momentum shell at 

— * 

| A;' | ~ A with the thickness t is given by 

id 



j p (jfe, w) = y^—S^ r d6 sin^ 2 6 ( 2D A 2 + —k 2 + 2DA U + — A^ 2 k 2 
{2-K) d T J Q \ 2 4 

+ ^(a-2)A a ~ 2 k 2 cos 2 9-2a-iTtu^ \ (A.17) 

where = 2i\ d l 2 jY[dj2) is the surface of <i-dimensional unit sphere. Carrying out the 
integration one obtains 

PA d K J ( 1 D + ^A CT - 2 + 122.(2- - l)A a ~ 2 ir \ 

'^-^r\h- — — k + ■ (A ' 18) 

where K d = Sd/(2ir) d and M = DA 2 + DA a — a. Thus, integrating out the one-loop 
corrections within the momentum shell one obtains to lowest order three terms which 
have the same structure as the terms in the free part of the action. Similarly we can 
compute the vertex loop integral, giving 

/c\ d kr\i,i 
^ 5fr Gg(fc,a;)Go(-fc,-a;) (A.19) 

_ J d d kdu 1 



i r d d i 



(2n) d + 1 (Dk 2 + D\k\° -a- iruj) 2 (Dk 2 + D\k\° -a- ituj) 
d d k 1 lK A A d 



> (27r) d A(Dk 2 + D\k\° - a) 2 ' AtM 2 



Renormalization group flow equations 

Absorbing the contributions from the shell integrals in the coefficients of the action, the 
renormalization group equations to one-loop order are given by 

de In r = - % - 2S 2 
di In a = z — r]^ — S\ 

a,,n D = ,-^-2- S2 (l-I(f + g-9SA- 

di In D = Z — T)ip — U (no loop corrections) (A. 20) 

d e hig = z - 3^/2 - d/2 - SS 2 

di In po = — (Vi> Vpo ) ' no loop corrccti ° ns ) 

where 

_ g 2 A d K d _ g 2 A d K d 

Si -^taT' S2 -J6tW (A - 21) 

As discussed in the main text, the coefficients D and po are not renormalized by loop 
corrections to any order of perturbation theory, implying the well-known generalized 
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hyperscaling relation for the critical initial slip exponent 

= Vp0 ~ z V ^ = - - 25 (A.22) 
as well as a new scaling relation among the three standard exponents 

-u ± (a-d) -2/3 = 0. (A.23) 
Fixed point and critical exponents 

Let us now assume that the coefficients D and r are constant during renormalization. 
Denoting by e = d c — d = 2a — d the difference from the upper critical dimension, the 
remaining RG flow equations take the form 

de In a = o — S\ 



cr-2 



d e lng = e/2-7S 2 
These equations have the fixed point 

S* = <r, S 2 * = e/4, (A.24) 

or in terms of the original coefficients 

2DA a 



7a 



-e + 0[t' 



(g*) 2 = 7 £ e + Otf) (A.25) 

In the vicinity of the fixed point the linearized RG flow is governed by the matrix 

, d a a{a-S x ) d a a{e/2-7D 2 ) . 
' d^a-SO d g g(e/2-7D 2 ) j 

' 9=9 ; a=a* 

To first order in e, the eigenvalues of this matrix are ai = a — 2e/7 and a 2 = — e. Here 
the first eigenvalue is positive and represents the attractive line of the renormalization 
group flow, meaning that a vanishes as A~ Ql as the fixed point is approached. Since 
a plays the role of the reduced percolation probability, we may identify u± with a^ 1 . 
Having reached the fixed point, the postulated stationarity of D and r implies that 
T]^ = — e/7 and z = a — e/7. Expressing these identities in terms of the standard expo- 
nents one is led to the results in Eq. (l4"3l) . 
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